Duality between equilibrium and nonequilibrium networks 
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In statistical physics any given system can be either at an equilibrium or away from it. Networks 
are not an exception. Most network models can be safely classified as either equilibrium or growing. 
Here we show that if certain conditions are satisfied, there exists a growing counterpart for any 
equilibrium network model, and vice versa. The equivalence between the two systems is exact not 
only asymptotically, but even for any finite system size. The required conditions are satisfied in 
causal sets and to a large extent in some real complex networks. 
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Statistical physics studies equilibrium and nonequi- 
librium systems using different theory and methods, as 
these systems are drastically different. Networks are not 
an exception. The vast realm of network models can 
be roughly divided into equilibrium and nonequilibrium 
domains [TJ [5] . In the former, one usually builds an equi- 
librium ensemble of graphs of a given size. Classic ex- 
amples are classical random graphs [3], soft configura- 
tion model and hidden variable models [4], or random 
geometric graphs [5]. In nonequilibrium models, graphs 
usually grow by adding nodes one at a time, introducing 
statistical dependencies. Preferential attachment [HI [7] is 
perhaps the best known example. These two approaches 
are clearly different [HIH11H]- ^ n the simplest example, 
i.e. classical random graphs Qn, p , each pair of N nodes 
are independently connected with the same probability 
p. The resulting degree distribution is Poissonian with 
mean k = pN. In a growing version of this model, new 
nodes N — 1, 2, 3, . . . are coming one at a time connect- 
ing to each existing nodes with probability p = k/N. The 
degree distribution is exponential pQ. 

Causal sets (or causets) [TO] can, on the one hand, 
be considered as random geometric graphs in Lorentzian 
spaces: distribute nodes uniformly at random, using 
a Poisson point process, over (a compact patch of) a 
Lorentzian manifold, and then connect each node pair 
located at timelike distance [10]. This construction is 
clearly equilibrium. On the other hand, causets were 
proposed as a model of quantum gravity, in which case 
they must be a dynamic model of growing spacetime [TT] . 
To approximate the spacetime that we observe, this cur- 
rently unknown dynamics, if run to completion (thermo- 
dynamic or large-size limit), must lead to asymptotically 
the same graphs as in the static case, reminiscent of the 
connection between Cayley trees and Bethe lattices [12] . 
This imperative duality between stasis and dynamics is 
somewhat unusual in physics. It raises many questions, 
including philosophical [13] . while the causet quantum 
gravity program is far from complete. 

Notwithstanding these difficulties, causets in asymp- 
totically de Sitter spacetimes [14], such as the spacetime 
of our accelerating universe [TU [16], are structurally 
similar to real complex networks, e.g. the brain or In- 



ternet [17j . This similarity is a consequence of asymp- 
totic equivalence between the growth dynamics of de Sit- 
ter causets and complex networks [17]. Specifically, a 
growing model of complex networks [18] is asymptotically 
identical to growing de Sitter causets [17], leading to an 
apparent paradox: given that only an equilibrium con- 
struction for de Sitter causets is currently known, what 
are these "growing causets," or more generally, how can 
a graph model be equilibrium and nonequilibrium at the 
same time? 

Here we show that if certain conditions are satisfied, 
then there exists a dynamic (growing) counterpart Q d for 
any static (equilibrium) graph model Qs, and vice versa. 
Specifically we prove that Qs and Qd are identical, and 
that this equivalence is exact not only asymptotically, 
but even for any finite graph size. That is, if Qs gener- 
ates graph G with probability P(G), then so does Qd- 
For concreteness we work with de Sitter causets in 1 + 1 
dimensions, and generalize at the end. 

Static 1 + 1-dimcnsional de Sitter causets can be de- 
fined as follows 17 . Given volume V and sprinkling den- 
sity 5 as parameters defining the expected size N = SV 
and average degree of resulting graphs: 

1) sample graph size N from the Poisson distribution 
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2) sample N random numbers Vi, i — 1,...,N, 
from the uniform distribution on [0, V], shorthand Vi £ 
U(0,V); 

3) assign to each node i its conformal time and space 
coordinates given by r/i — arctan[wi/(27r)] and 0i € 
W(0,2tt); 

4) connect each pair of nodes i and j if Ai] > AO, 
where A 77 = \rji — r)j | and AO = it — \tt — \0i — 3 \ | are the 
temporal and spatial distances between the two nodes. 

The first two steps in this definition are nothing but 
an implementation of a Poisson point process (PPP) with 
rate 5 on [0, V"]. The third step is a PPP on a compact 
patch V of volume V in de Sitter spacetime between con- 
formal times r\ = and 770 = arctan[V/(27r)]. The ele- 
ment of volume in this spacetime is dV = sec 2 rj drj dO, 
and step (3) ensures that the spatial and temporal node 
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FIG. 1: Illustration of the Poisson point process in Gd,v (top) 
and Qd.n (bottom). 



densities are p(9) = l/(2ir) and p{rj) = sec 2 r\j tan r)o, 
meaning that nodes are indeed distribution uniformly at 
random in V according to the volume form dV. The re- 
sulting graphs have a power-law distribution of node de- 
grees k, P(k) ~ fc~ 7 , with 7 = 2, and strongest possible 
clustering [17 . We denote this ensemble by Gsy- 

We now want to build a dynamic version of this ensem- 
ble, Gd,v, by adding nodes one at a time. Unfortunately 
it is impossible if we want Gd,v to be identical to Gsy- 
The best we can do is the following Gd v definition. For 
7 = 1,2,...: 

1) select any real number Vj > and let V = J2j=i Vj; 

2) sample iV/ € P«vj(iV» and let N = Y?j=i N J> 

3) for i = N-Nj + 1,...,N, sample Vi £ U{V-V U V) 
and 9i € U(0, 2ir), and set r]i = arctan[v i /(27r)]; 

4) connect all new nodes i = N — Ni + 1, . . . , N to all 
existing nodes j = 1, . . . , N if A 77 > A#. 

Each time / we simply grow the spacetime patch V that 
nodes occupy, by adding to V a new spacetime chunk 
Vi of volume Vj, and throwing Nj new nodes into Vi, 
Fig. [TJtop). If all Vi — 1/8, then the expected num- 
ber of nodes we add each time is Nj = 1. However, we 
cannot force the number of new nodes Nj to be exactly 
1. Indeed, suppose we freeze the Gd,v growing process 
at some volume V , and compare the result to the static 
Gs,v construction. If we have been adding exactly one 
node in each V/, then the number of nodes Ni = 1 in 
Vi is clearly different from what it is in Gsy, where this 
number can be anything — by the PPP definition, it is 
Poisson-distributed around 1. Step (2) in the Gd,v defi- 
nition ensures that the JVj distributions are the same in 
both cases. The expected numbers of nodes N in gener- 
ated graphs G are also the same, for any partition of V in 
Gsy into chunks {V/} in Gd.v, because a sum of Poisson- 
distributed random numbers is Poisson-distributed with 
the mean equal to the sum of means. 

Probability P(G) to generate any graph G is the same 
m Gsy an d Gd,v, i-e. the two ensembles are identical, if 
the linking conditions for node pairs are the same, which 
is clearly the case here (Ary > A9), and if the distribu- 
tions of node coordinates (77, 9) are exactly the same, too. 
The latter does hold as a simple consequence of the PPP 
definition, i.e. of the fact the expected number of nodes 
in any volume element dV is dN — 8 dV. For example, 
let N be the number of nodes in a Gsy sample. Then 
in the u-coordinates, the expected number of nodes lying 
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FIG. 2: Probability density functions (PDFs) for conformal 
time coordinates r\ in static and dynamic ensembles. Pan- 
els (a) and (b) deal with the fixed- V ensembles Gsy and 
Qd,v, where the exact 77-PDF is given by Eq. jSJ, the solid 
blue curves. The red circles and black dots are simulation 
results for Gsy and Gd,v- In (a), V = 2 x 2tt, 1 = 2, 
Vi = V 2 = 2n, and 770 = 1.107148. In (b), V = 2vr x 10 4 , 
I = 10 4 , Vi = . . . = V W i = 2tt, and 770 = 1.570696. Panels (c) 
and (d) show the corresponding results for fixed- iV ensembles 
Gs,n and Gd.n, where the exact PDF is given by Eq. d9J. 
In (c), AT = 2, I = 2, Ni = N 2 = 1. In (d), N = 10 4 , 
/ = 10 4 , Ni = ... = N 10 4 = 1. The numbers of samples 
(PPP runs) in simulations are S = 10 6 in (a) and (c), and 
S = 10 4 in (b) and (d). The sprinkling density (PPP rate) is 
8=1 everywhere. 



between v and v + dv in the sample is Np(v) dv, where 
p(v) = 1/V is the PDF of w-coordinates. The expected 
number of nodes in [v,v + dv] in the static ensemble is 
then 



dN s {v) = V Np(v)V S v(N)dv = 8dv. 



N=0 



(2) 



Assume now v £ Vi for some I in Gd,v, where the v- 
PDF is pi(v) — 1/V/. The same number in the dynamic 
ensemble is then 

00 

dN D (v) = Njpj^PsvjiNj) dv = Sdv, (3) 

Ni=0 

i.e. dN D (v) = dN s (v). As a consequence, in the rj- 
coordinates, the expected number of nodes in an infinites- 
imal spacetime chunk between rj and 77 + drj is in both 
cases 

dN(r)) = 8 dV(<n) = 2nS sec 2 r\ drj, (4) 

so that the probability density functions (PDFs) for 77 
are indeed the same, Fig. [2]ja,b), and equal to 

dN(rj) 



sec 2 rj 



N dr\ tan 770 



(5) 
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The equivalence between the coordinate one-point 
PDFs p(rj) is a necessary but not sufficient condition for 
two ensembles to be the same. The sufficient condition 
is that the joint distributions p(r}i,r]2, . . .) are also the 
same. If they are different, then even if their marginals 
p(r)) are the same, some other graph properties can be 
different, even in the thermodynamic limit. One example 
would be the order statistics jTHj defining, among other 
things, the distribution of the smallest time coordinate 
?7min of nodes, and consequently the distributions of their 
largest fc max and average k degrees. The joint coordinate 
distributions in Gs,v and Gd.v are the same, also as a 
PPP consequence. Both ensembles implement the same 
PPP process, so that due to statistical independence of 
PPP events (nodes), the joint PDFs are products of one- 
point PDFs in both cases. In other words, the two en- 
sembles are indeed exactly the same, and generate any 
graph G with the same probability P(G). 

But is it really impossible to add exactly one, or to that 
end, any fixed number of nodes at a time, as in many pop- 
ular models of growing networks, including preferential 
attachment? The answer is that it is certainly possible, 
but the resulting ensembles Gs,n and Gd.n that we de- 
fine next are identical to Gs,v and Gd,v only in the ther- 
modynamic limit. Yet in a certain sense that becomes 
evident from the exposition below, Gs,n and Gd,n are 
dual to Gs,v and Gd.v for any finite graph size. 

We first define the static ensemble Gs,n of graphs of 
a fixed size N, dual to Gs,v- Given graph size N and 
sprinkling density S as parameters defining the expected 
spacetime volume V = N/S occupied by nodes, and the 
average degree in resulting graphs: 

1) sample spacetime volume V from the Gamma dis- 
tribution 
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dual to the Poisson distribution, and set vn — V; 

2) sample N — 1 random numbers u, £ U(0,V), i = 
1,...,N-1; 

3) assign to each node i = 1, . . . , N its conformal time 
and space coordinates given by rji — arct an [7^/(271")] and 
i eU(Q/2n): 

4) connect each pair of nodes i and j if A77 > A9. 
Since the Gamma distribution T^ t g(V) is the distribution 
of waiting time V for the JV'th outcome of a PPP with 
rate 6, the Gs,n definition implements the same PPP as 
Gs,v , except that not the waiting time V but the number 
of events N is now fixed. 

The dynamic variant 5dm is now obvious. For / = 
1,2,...: 

1) select any integer number Nj > and let N = 

2) sample Vj £ Tn Iv s(Vi), let V — Ylj=i Vj> an d set 
v N = V; 

3) for i = N - Nj + 1, . . . , N - 1, sample v, £ 
U(V - Vj, V), and for i = N - N r + 1, . . . ,N sample 
9i £ U(0, 2ir) and set r\i = arctan[w i /(27r)]; 



4) connect all new nodes i = N — Ni + 1, . . . , N to all 
existing nodes j = 1, . . . , N if Ary > A9. 
As in the fixed- V case, we add nodes in bunches, except 
that the number of new nodes that we add each time is 
now fixed, see Fig.JIJbottom). In particular this number 
can be always 1. Yet we cannot strictly control the vol- 
ume that these nodes occupy, as we could not control the 
number of nodes in the volume-controlled case G*,v- 

The Gs,n and Gd.n ensembles are identical for any 
partition {Nj} of any N, since the two ensembles are 
two different implementations of the same fixed- TV PPP. 
In particular, the expected volume V that nodes occupy 
is the same in both cases because the sum of Gamma- 
distributed random variables is Gamma-distributed with 
the mean equal to the sum of means. The expected num- 
ber of nodes in any volume element dV is dN = 8 dV as 
before, but the one-point PDFs for time coordinates are 
different from G*y, where the v- and ?7-PDFs are defined 
on finite intervals [0, V] and [0,770] with 770 < tt/2, while 
in G*,n even at N = 1 the single node can have an ar- 
bitrarily large coordinate, meaning that the correspond- 
ing distributions are defined on the whole infinite space 
[0,oo) and [0, 7r/2). For example, since the distribution 
of the w-coordinate of i'th node in the PPP is given by 
Pi(v) — Tij(v), the PDF of ^-coordinates for N nodes is 
simply 



p{v) 
Q(N,Sv) 



1 N 

i=l 

T(N, 5v) 

r(A0 
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is the regularized Gamma function. The 77-PDF is then 



p(v) 



N 



sec 2 77 Q(N, 5 tan 77), 



(9) 



see Fig. §c,d). 

It is instructive to consider the following ensemble 
Gw,n where 'W stands for "wrong": for i — 1,...,N, 
sample Vi £ Ti g(Vi), and then finish the graph construc- 
tion as before. The intuition might be that since i'th 
coordinate is distributed according to r^^u), we can 
safely sample directly from this distribution, and one 
can check that if we do so, then the 7J-PDF is as needed, 
p(v) = (6/N)Q(N,Sv). However, this reasoning is wrong. 
The process is no longer a PPP, and the joint PDFs are 
different from Gs.n and Gd,n, because the ordering of 
coordinates can now be violated, e.g. Vi can with cer- 
tain probability be larger than Vj for any i < j, while by 
definition, Vi t g(v) is the distribution of the i'th largest 
coordinate in a PPP. As a result Gw,n is not identical to 
Gs,n and Gd.n even in the thermodynamic limit, which 
one can verify in simulations by looking at the maximum 
and average degree statistics, for instance. 

In the thermodynamic limit N — > 00, the regularized 
Gamma function approaches I, Q(N,Sv) —> 1, so that 
p(v) — > S/N and ^(77) — > (S/N) sec 2 77, and the fixed-A^ 
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FIG. 3: Degree distribution (a) and clustering (b) in de 
Sitter causets Gs,n and Gd.n, 1000 samples each of size TV = 
23752, sprinkling density 8 = 0. 83475 1/2-7T, average degree 
k = 4.9188 and average clustering c = 0.7862 and c = 0.7854, 
respectively, and in some real networks: AS Internet (TV = 
23752, k = 4.9188, c = 0.6055), functional brain network 
(N = 23713, k = 6.1436, c = 0.1552), and PGP Web of Trust 
(TV = 23797, k = 7.8587, c = 0.4816). 

ensembles Gs,n a* n( I Gd,n become asymptotically identi- 
cal to the fixed- V ensembles Gs,V and Gd,v, a U imple- 
menting the same PPP on the infinite upper half of de 
Sitter spacetime with 77 > 0. In other words, we have 
proved the following diagram 

Gs.y = Gd,v 

I I 
Gs,n = Gd,n 

where the equal sign '=' means the exact equivalence 
for any system size, while symbol £ ~' stands for asymp- 
totic equivalence and for the V.vs.iV PPP duality at finite 
sizes. 

In general, the same picture would apply to any pair of 
static and dynamic network models in which the joint dis- 
tributions of (hidden or latent) variables |4 , attributes, 
or coordinates of nodes, including their birth times, are 
the same, and in which the probabilities of connections 
between nodes, as functions of these hidden variables, 
are exactly the same, too. Surprisingly, these two condi- 
tions are rarely satisfied. As a result, static and evolving 
networks are analyzed using quite different methods. To 
the best of our knowledge, there is no equilibrium con- 



struction that would be exactly identical to preferential 
attachment, for instance. Yet the required conditions are 
satisfied in causal sets and more generally, in random 
geometric graphs, including random geometric graphs 
in hyperbolic spaces |18j . These hyperbolic graphs are 
asymptotically identical to de Sitter causets [17] , and de- 
scribe well not only the structure but also the dynam- 
ics of some real networks [18]. Here we have to stress 
that this hyperbolic-de Sitter equivalence is exact only 
for a specific (default) set of parameters [18] . leading to 
specific types of graphs — namely, to graphs with power- 
law exponent 7 = 2, and with strongest possible clus- 
tering, i.e. zero temperature [18] . The former property 
applies to many real networks, Fig. [3^a), but the latter 
does not, explaining differences between clustering in de 
Sitter causets and real networks, Fig.[3jb). However, the 
temperature of real networks is low, meaning that their 
clustering, while not the strongest possible, is still strong 
and thermodynamically stable [T8J. In that respect the 
differences between de Sitter causets and real networks 
may be not that big. 

In conclusion, almost all real networks are growing, 
justifying certain concerns, if not skepticism, about the 
utility of equilibrium methods in analyzing real networks. 
Yet the structure and dynamics of these networks turn 
out to be well described by network models characterized 
by the unusual exact equivalence between equilibrium 
and nonequilibrium formulations that we have proved 
here. These results thus provide a different perspective 
and further theoretical grounds for the use of power- 
ful equilibrium methods, including exponential random 
graphs and other graph entropy tools 20 2!j] , in the anal- 
ysis of real networks. 
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